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ABSTRACT 

Context. The highly dynamic atmosphere above sunspots exhibits a wealth of magnetohydrodynamic (MHD) waves. Recent studies 
suggest a coupled nature of the most prominent phenomena: umbral flashes and running penumbral waves (RPWs). 

Aims. From an observational point of view, we perform a height-dependent study of RPWs, compare their wave characteristics, and 
aim to track down these so far only chromospherically observed phenomena to photospheric layers to prove the upward propagating 
field-guided nature of RPWs. 

Methods. We analyze a time series (58 min) of multiwavelength observations of an isolated circular sunspot (NOAA11823) taken at 
high spatial and temporal resolution in spectroscopic mode with the Interferometric Bldimensional Spectro-polarimeter (IBIS/DST). 
By means of a multilayer intensity sampling, velocity comparisons, wavelet power analysis, and sectorial studies of time slices, we 
retrieve the power distribution, characteristic periodicities, and propagation characteristics of sunspot waves at photospheric and chro¬ 
mospheric levels. 

Results. Signatures of RPWs are found at photospheric layers. Those continuous oscillations occur preferably at periods between 
4-6 min starting at the inner penumbral boundary. The photospheric oscillations all have a slightly delayed, more defined chromo¬ 
spheric counterpart with larger relative velocities, which are linked to preceding umbral flash events. In all of the layers, the power 
of RPWs follows a filamentary fine-structure and shows a typical ring-shaped power distribution increasing in radius for larger wave 
periods. The analysis of time slices reveals apparent horizontal velocities for RPWs at photospheric layers of « 51 km/s, which de¬ 
crease to « 37 km/s at chromospheric heights. The photospheric distribution of peak periods at the inner penumbra resembles the 
chromospheric cases. 

Conclusions. The observations strongly support the scenario of RPWs being upward propagating slow-mode waves guided by the 
magnetic field lines. Clear evidence for RPWs at photospheric layers is given. Assuming an inverse proportionality of the peak pe¬ 
riod and cut-off period on the cosine of the field inclination, the penumbral magnetic field inclination is increasing toward the outer 
penumbra. The more rapid increase and the larger horizontal velocities at photospheric heights hint at the more horizontal penumbral 
field inclination at lower heights. 
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1. Introduction 

Sunspot waves are one of the most spectacular and dynamic phe¬ 
nomena in the solar atmosphere. Since their first detection, the 
question about the nature of umbral flashe s ([Beckers & Tallant 
1969 ) and r unning penumbral waves ( Giovanelli|1972| Zirin & 
Stein| | 19721 henceforth RPWs) and their relation has been un¬ 


der intense discussion. Whereas umbral flashes have been com¬ 
monly interpreted as upward propagating magnetoacoustic slow¬ 
mode waves guided by magnetic field lin es (e.g., |Centeno et al.| 
|2006[ |de la Cruz Rodriguez et al.||2013| ), the case for RPWs is 
not yet clear. 

A first scenario suggests that RPWs are purely chromo¬ 
spheric waves excited by umbral flashes and propagate radially 
outwa rd across the chromospheric penumbra (Tziotz iou et al.| 
|2006| . Secondly, recent studies support the scenario that run¬ 
ning penumbral waves, like umbral flashes, are channeled, up¬ 


ward propagating waves excited at lower layers (Rouppe van 


der Voort et al. 

2003; Kobanov et al. 2006; Bloomfield et al. 

2007; Jess et al. 

2013 e.g.,). The increasing magnetic field in- 


clination toward the outer penumbra would explain this visual 
pattern of horizontal propagation at a chromospheric layer since 
the increasing path length for a coherent wave front delays its oc¬ 
currence. However, a detailed photospheric observation focusing 


on RPWs to disprove the purely chromospheric nature of RPWs 
was still lacking until now. 

Under the assumption of field-guided waves, the peak period 
or frequency in a power spectrum of sunspot waves can be used 
as an indicator for the magnetic field inclination. The acoustic 
cut-off, which reduces waves with periods above a certain cut¬ 
off value, is highly dependent on the field topology ( Bel & Leroy 


1977 ). The effective cut-off period increases with the inclination 
angle of the field lines. Observational studies (e.g., |Jess et al" ~ 


201 3] [Wan et al. ||20 1 4| ) show that the peak period in chromo¬ 


spheric power spectra increases from 3 min in the umbral area to 
around 10 min toward the outer penumbra. For further informa¬ 
tion, see Bogdan & Judge (2006). 

Most studies on RPWs have focused on the chromosphere. 
Only a few attempts (e.g., |Lites et al.|T998[ |Bloomfield et al. 
|2007| ) have been made to determine their photospheric behavior 
in detail. In this work, we perform a novel high-resolution, mul¬ 
tilayer study concentrating on the photospheric characteristics 
of RPWs. In Sect. [2] of this article, we give a brief description of 
the observations, the data, and calibration techniques. In Sect.[3j 
we present the wave analysis and discuss the results. Finally, in 
Sect.[4]we draw our conclusions and describe the importance of 
the new findings. 
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Fig. 1 . NOAA11823 as seen in continuum (left panels) and line-of-sight magnetic field strength (right panels) on August 21st 2013 
at 15:00:45UTC as observed by HMI (a+c+d) and ROSA (b). The black squares mark the analyzed region. The insets (b+d) show 
the sunspot with contours indicating the inner and outer penumbral boundary from the time-averaged HMI continuum intensity. The 
HMI magnetogram is scaled from -1 kG (black) to 2kG (white). 


2. Observations and data analysis 

To study the wave phenomena in a sunspot (see Fig.[l]) from 
photospheric to chromospheric layers, observations were carried 
out on August 21st 2013 from 14:53UTC to 15:51UTC at the 
Dunn Solar Telescope (DST) at the National Solar Observatory 
in New Mexico. The atmospheric conditions, which are crucial 
for ground-based observations, were excellent and stable allow¬ 
ing an effective spatial resolution of up to 0.4 " even for the non- 
reconstructed data. Simultaneous multi-instrument, multiwave¬ 
length observations were performed with the etalon-based imag¬ 
ing spectropolarimeter IBIS ( |Cavallini]|2006| ) in spectrometric 
service mode toget her with the broadband imaging instrument 
ROSA (Jess et al. 2010[ ) to sample the sunspot’s photosphere 
and chromosphere (see left panels of Fig. [2]). A complete and de¬ 
tailed overview of the data will be given in a forthcoming article, 
including a full time lapse movie of the observations (Lohner- 
Bottcher & Bello Gonzalez, in preparation). 

In this work, we focus on the narrowband spectromet¬ 
ric observations done with the Interferometric Bldimensional 
Spectro-polarimeter (IBIS) for the spectral lines Fel 630.15 nm, 
Na ID1 589.6 nm and Ca II 854.2 nm. For the penumbra, the Fe I 
630.15 nm intensities stem from photospheric height of 0 km 
(line continuum) to 300 km (line minimum) above Tsoonm = 1 
( [Bello Gonzalez et al.|2005| ). Beyond that, the NalDl 589.6nm 
and Call 854.2nm lines cover the sunspot’s atmosphere up to 
chromospheric layers. Each spectral line was sampled at 10 - 12 
nonequidistant wavelength positions (exposure time of 35 ms). 
The lines were scanned successively, leading to a cadence of 
13.2 s. According to the Nyquist criterion, the oscillations with 
periods down to 26 s can be investigated. 

The isolated, perfectly stable sunspot of NOAA11823 is 
shown in Fig. [T] in continuum intensity (panels a+b) and line- 
of-sight magnetic field strength (panels c+d). The sunspot has 
a fully-developed circular penumbra, a diameter of 24 Mm and 
is located close to disk center at a heliocentric angle 6 = 14° 
at disk coordinates (X, Y) = (63",-222"). The HMI magne¬ 
tograms revealed the unipolar direction of the magnetic field and 
a maximum umbral field strength of +2.2 kG. As continuum in¬ 


tensity and magnetic field strength show, some small and stable 
pores surround the sunspot at a distance of more than 30" from 
the outer penumbral boundary. The reduced field of view, indi¬ 
cated by Figs.[l]b+d, has a size of 42 Mm 2 at a pixel scale of 
0.098 "px” 1 . 

The data calibration of the IBIS data involved background 
intensity subtraction, flat-field calibration, correction for colli¬ 
mated wavelength shifts and prefilter transmission, as well as 
an iterative reduction of temporal image distortions caused by 
atmospheric turbulences. Doppler velocities for Fel 630.15nm 
and NalDl 589.6nm were calculated as relative shifts of line 
minimum positions of a Gaussian approximation to the full line 
profile. The calculated velocities therefore represent the aver¬ 
aged motion over the atmospheric layers these lines form in. 
For the penumbra, this yields reliable results since the oscil¬ 
latory line shifts dominate the effects of the magnetic field on 
the line profiles. The accuracy of the photospheric velocities 
is confirmed by the Doppler velocities from the Helioseismic 
and Magnetic Imager (HMI/SDO). These continuous and stable 
satellite data using the Fel line at 617.3nm with a cadence of 
45 s and spatial scale of 0.5 "px -1 serve as context information 
for our observations (see Fig.§. 

For the investigation of characteristic frequencies and power 
distribution, we performed a time-position-dependent wave 
powe r analys is based on wavelet techniq ues ( |Torrence~& Compo 
|1998| ) as in |Bello Gonzalez et al.| ( |201Q] ). On the basis of a con¬ 
tinuous oscillatory stability, we limit this study to the interpre¬ 
tation of the global (time-averaged) wavelet power spectra cal¬ 
culated at a fine nonlinear scale sampling. One hundred period 
steps sample the range from 26 s to 14 min; the step size in¬ 
creases to higher periods. The characteristic interval of sunspot 
waves between 2 min and 6 min is scanned by 34 steps. For the 
power analysis with respect to the distance to the sunspot center, 
we included spherical projection effects implied by the sunspot 
location and orientation. We therefore used foreshortened cir¬ 
cles to compute the azimuthal average as in Lohner-Bottcher & 
Schlichenmaier] ( |2013| ). Close to disk center, these effects are 
small but should not be neglected. For the temporal evolution of 
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Fig. 2. Three-dimensional view of intensities and peak periods of intensity wave power of NOAA11823 at various wavelength 
positions. The intensities (left) show the sunspot at August 21st 2013 at 15:00:06UTC (±3 s). The images along the z-axis belong to 
several line core and wing positions of Fel 630.15 nm, Nal 589.6nm, and Call 854.2nm. The corresponding time-averaged (« lh) 
distribution of peak periods Tpeak of the intensity wave power is shown on the right. The periods are scaled from 2.5 min (dark 
blue) to 8 min (dark red). The black contours indicate the location of the umbra (inner) and penumbra (outer) in continuum intensity 
(bottom panel). Whereas the length of the axis arrows represent distances around 1.5 Mm, the image positions along the z-axis are 
not to scale. 


RPW velocities in radial direction, we also used the projection 
of circular sectors. 


3. Results and discussion 

The analysis of this high-resolution, multiwavelength observa¬ 
tion of a perfectly stable and symmetrical sunspot yields clear 
evidence for photospheric RPWs and their field-guided propa¬ 
gation to higher layers. 

The wave power analysis ( P ) for Doppler velocities from the 
photospheric Fel 630.15nm and photospheric-chromospheric 
NalDl 589.6 nm lines reveals the spatial distribution and mag¬ 
nitude of the oscillatory features in the sunspot’s atmosphere. 
As shown in Fig.[3| the power of the RPWs is concentrated in 
the inner penumbra for both layers (bright colors indicate high 
power). Whereas the first RPW signatures are detected at the 
umbral boundary at periods from 3-4 min, the main power peaks 
between 4-6 min as filamentary structures in the inner half of 
the penumbra. The averaged power in this period range is shown 
in Fig.[3]a and b. The mentioned fine-structured pattern can also 
be seen in the calculated peak period maps (Figs. [2] and [4]), and 
seems to reflect the penumbral spine-intraspine distribution in 
intensity. 

As the illustrative contours in Fig.[3]c and d point out, with 
increasing distance in the inner penumbra the RPW power Px 
shifts toward larger periods. The shift in period at the inner 


penumbra stands in contrast to the constant /7-mode signal in 
the 4-6 min range at the outer penumbra and quiet sun. This 
evidence for RPWs is even more prominent in the low chro¬ 
mosphere where the shift to higher periods is still traceable to 
the outer penumbra up to periods of around 10 min. To verify 
this radial increase of wave periods, which for chromospheric 
heights is in line with recent studies ( Jes s et al.|2013|fYfian et al. 
2014), the peak periods in the power spectra of both velocity 
(Fig.|4]) and intensity (Fig. [2]) oscillations were determined for 
all heights. Whereas at photospheric layers the radial transition 
to larger peak periods (from 4 min to 8 min) is centered in the 
innermost penumbra, at higher layers the increase in period is 
smoother. In chromospheric heights, larger periods are therefore 
located further outside in the penumbra and beyond. 


The height-sampled distribution of wave power and shift 
in peak period hint at the opening field topology of s unspots. 
According to cut-off theory and observations (|Bel & Leroy11977 [ 


De Pontieu et al.|2004[|Tziotziou et al.|2006| ), the observed peak 

period r Pea k of field-guided waves depends on the acoustic cut¬ 
off period T cut ^ and increases with temperature #k and incli¬ 
nation angle 0 (with respect to the vertical) qualitatively like 
Tpeak ~ T cut , 0 /1.25 ~ V^k / cos 0-Larger inclinations would 
especially allow the propagation of overpowering high-period 
waves into higher atmospheric layers. Assuming the acoustic 
cut-off already occurring for photospheric RPWs, the larger pe¬ 
riods would indicate a more horizontal orientation of the penum- 
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Fig. 3. Spatial distribution of the time-averaged wavelet power P 
of the Doppler velocities from Nal 589.6 nm (upper panels) and 
Fe I 630.15 nm (lower panels). Panels a and b show the power in 
the 4-6 min band. The power (bright colors indicate high power) 
is plotted in logarithmic (arbitrary) units to assure best contrast 
in the penumbra. The contours mark the inner and outer penum¬ 
bral boundary in continuum intensity. Circular sectors 1-3 (white 
solid) follow the foreshortened projection of a circle (dotted). 
Panels c and d show the azimuthally averaged global power spec¬ 
trum P x for periods T (in min) according to the distance X (in 
arcsec) from the spot center. The white contours mark the impor¬ 
tant elements in the spatial power analysis. For spatial compari¬ 
son, an arbitrary radial element in continuum intensity is added 
below. The dashed lines indicate the inner and outer penumbral 
boundaries. 


bra at photospheric than at chromospheric levels. In addition 
to the spectroscopic study of IBIS data, the analysis of HMI 
Doppler velocities and their power distribution confirm the pho¬ 
tospheric signatures of RPWs, though at a lower spatial and tem¬ 
poral resolution (compare Fig.[4]a+b). 

To prove the field-aligned propagation of RPWs in the lower 
sunspot atmosphere, we study the temporal evolution and trajec¬ 
tories of the waves. Therefore, we perform a sectorial time-slice 
analysis for the photospheric and chromospheric Doppler veloc¬ 
ities, as shown in Fig. [5] To investigate the fluctuations, we sub¬ 
tract the temporally averaged Doppler velocities and focus on the 
relative velocities. Since the RPWs seem to follow the filamen¬ 
tary fine-structure of the penumbra and under the assumption 
of coherent wave trains for nearby fields, we select three circu¬ 
lar (foreshortened) sectors with opening angles of 10° for the 
penumbral regions with highest velocity power (see Fig.[3]a,b). 
The temporal behavior of the azimuthal average at each distance 
from the spot center is shown in Fig. [5] for the velocities derived 
from Nal 589.6nm (panels a-c) and Fel 630.15nm (panels d- 
f). A detailed study of the oscillatory pattern demonstrates the 
photospheric origin of RPWs as follows: 

- Continuous and clear signatures of penumbral oscillation at 
photospheric height are found. 

- Photospheric oscillations at the inner penumbra have an am¬ 
plitude of up to 0.3 km/s. 


05 

O 

CO 


_05 

O 

CO 


-240 

-230 

-220 

-210 

-240 

-230 

-220 

-210 



50 60 70 80 

Solar-X ["] 


Fig. 4. Spatial distribution of the peak periods from the global 
power spectra of Doppler velocity oscillations. The major pe¬ 
riods Tpeak (in min) are shown for a) Fel 630.15nm, b) Fel 
617.3nm (HMI), and c) Nal 589.6nm. The scale ranges from 
2.5 min (dark blue) to 8 min (dark red). The black contours mark 
the umbral and penumbral boundaries in continuum intensity. 
The white arrow is pointing to the disk center. 


- All sectors provide observational evidence for outward- 
directed horizontal propagation in the inner penumbra. 

- In the inner penumbra close to the umbral boundary, all 
photospheric waves have a slightly delayed chromospheric 
counterpart with larger oscillatory amplitudes of up to 
0.8 km/s. 

- A simple analysis of the slopes of all wave trains over the 
first 4" of the inner penumbra was performed in the sectorial 
time slices shown in Fig. [5] The analysis reveals the apparent 
horizontal velocities, Vror, combined to Fig. [6] For the pho¬ 
tospheric layers (black symbols and histogram), represented 
by the Fe I 630.15 nm line, the apparent horizontal velocities 
average to (vhor,fci) = 51 ± 13 km/s. 

- At upper photospheric to chromospheric height, represented 
by the Nal 589.6nm line (Fig.|5j red symbols and his¬ 
togram), the wave trains indicate smaller apparent horizontal 
velocities of <VHOR,Nai) = 37 ± 10 km/s. The smaller veloc¬ 
ity is illustrated exemplarily in Fig. [5] by the steeper slopes 
compared to the photospheric case. 


At chromospheric layers, the retrieved apparent horizontal 
velocities of around 20 - 50 km/s at the inner penumbra are in 
line with recent studies (e.g. , [Tziotziou et al.|2006 Kobanov 


et al.||2006} |Jess et alT]|2013| ). Commonly, this apparent propa¬ 


gation speed decreases radially toward the outer penumbra. As 
observations at higher chromospheric to transition region layers 
have shown ( [Madsen et al.|2015) , the retrieved horizontal veloc¬ 
ities decrease to around 10 km/s. Vice versa, the larger appar¬ 
ent velocities of around 30 - 70 km/s at photospheric layers fit 
well into the trend of decreasing horizontal velocities at the in¬ 
ner penumbra toward higher layers. The topological model of the 
sunspot’s magnetic field (e.g., jWestendorp Plaza et al. 1997 ), the 


4 






















































J. Lohner-Bottcher and N. Bello Gonzalez: Signatures of running penumbral waves in sunspot photospheres 


Sector 1 Sector 2 Sector 3 



X ["] X ["] X [ n ; 


Fig. 5. Temporal evolution (in min) of the relative velocities in 
the sunspot atmosphere. The analysis of time slices was per¬ 
formed for the Doppler velocities VNai (panels a-c) of Nal 
589.6 nm and Vp e i (d-f) of Fel 630.15 nm in all three sectors, 
also shown in Fig. [3] For each sector, the azimuthal average at a 
radial distance X (in arcsec) from the spot center was calculated. 
The dashed lines mark the umbral boundary. The velocity scale 
AV LO s ranges between ±0.6 km/s (a-c) and ±0.3 km/s (d-f). The 
black bars trace the apparent wave trains. 


field-guided propagation of running penumbral waves, and their 
visual appearance at a certain layer can explain this behavior 
( [Bloomfield et al.[2007 [[Madsen et al.|2015 ). In the penumbra, 
the magnetic field inclination increases radially from the ver¬ 
tical umbral field. At photospheric layers, this bending of the 
magnetic field lines is stronger than at higher atmospheric lay¬ 
ers in which the radial increase in field inclination happens more 
smoothly. As stated by Bogdan & Judge ( 2006), the large appar¬ 
ent horizontal velocity of RPWs (here up to 80 km/s in Fig. [6]) 
rather reflects relative travel time differences for waves guided 
by the individual less and more inclined magnetic field lines of 
the inner penumbra. 


4. Conclusions 


We provide strong evidence for the presence of RPWs in the 
sunspot photosphere. Our observations conflict with the scenario 
that RPWs are a purely chromospheric wave phenomenon and 
strongly support the recent theory (Rouppe van der Voort et al. 
2003[ |Bloomfield et al.|2007[|Jesse t al. 2013, e.g.,) that RPWs 
are upward propagating slow-mode waves guided by the mag¬ 
netic field lines. Under the assumption of an expanding magnetic 
field topology with more inclined field lines with radial distance 
from the vertically oriented umbra, the performed wave power 
analysis substantiates the dependence of the power spectra on 
the atmospheric height and inclination of the magnetic field. 



time t [min] #N 


Fig. 6. Apparent horizontal velocities, Vror, of running penum¬ 
bral waves at two atmospheric layers. The results for the photo- 
spheric Fel line (Fig.[5]d-f) are shown in black, for the chromo¬ 
spheric Nal line (Fig.]5]a-c) in red. Left panel: the apparent hor¬ 
izontal velocities (in km/s) over the first 4" of the inner penum¬ 
bra are plotted according to the observation time, t. The aster¬ 
isks, diamonds, and triangles show the apparent speeds of the 
wave trains in sectors 1, 2, and 3. The dashed lines indicate the 
average apparent velocities, (vhor,fci) and <VHOR,Nai) for both 
layers. The averages and their standard deviations are given in 
the figure legend. Right panel: histogram showing the number 
#N of velocity values vror from the left panel within a 5 km/s 
binning interval. 


The solar atmosphere exhibits a wealth of dynamical phe¬ 
nomena at all scales. These processes and their effects are 
strongly coupled. Especially in the case of sunspots, oscilla¬ 
tions, flows, and other dynamical effects can interact and have 


to be taken into account in observational studies (e.g., Esteban 
|Pozuelo et al.|2015] ). 

A forthcoming, extensive study of the presented data 
(Lohner-Bottcher & Bello Gonzalez, in preparation) will include 
an analysis of umbral flashes in the chromosphere (a bright um¬ 
bral flash event can be seen in Fig.[2]), their characteristic simi¬ 
larities and differences from RPWs, and the possible reconstruc¬ 
tion of the magnetic field topology using the wave characteris¬ 
tics in sunspots. To verify the existence of photospheric RPWs, 
we suggest further observations with high spatial and temporal 
resolution using photospheric spectral lines that are magnetic in- 
sensitive (with a Lande-facto r g=0), for examp le Fel 557.6nm 
( Bello Gonzalez et al.||2009| or Fel 543.4nm ( [Bello Gonzalez] 
et al.|2010| ). For the latter, we have also found RPW signatures. 


As shown in Fig. [4] signatures of running penumbral waves can 
also be found in HMI Dopplergrams. For a quantitative evalua¬ 
tion of sunspot waves, we propose a statistical evaluation using 
a large sample of sunspots. 
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